Sub-Poissonian statistics of Rydberg-interacting dark-state polaritons 
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Interfacing light and matter at the quantum level is at the heart of modern atomic and optical 
physics and enables new quantum technologies involving the manipulation of single photons and 
atoms. A prototypical atom-light interface is electromagnetically induced transparency, in which 
quantum interference gives rise to hybrid states of photons and atoms called dark- state polaritons. 
We have observed individual dark-state polaritons as they propagate through an ultracold atomic 
gas involving Rydberg states. Strong long-range interactions between Rydberg atoms give rise to an 
effective interaction blockade for dark-state polaritons, which results in large optical nonlinearities 
and modified polariton number statistics. The observed statistical fluctuations drop well below the 
quantum noise limit indicating that photon correlations modified by the strong interactions have a 
significant back-action on the Rydberg atom statistics. 



Electromagnetically induced transparency (EIT) is a 
quantum interference effect in which an otherwise absorb- 
ing medium is rendered transparent for a weak probe field 
by an optical control field [1 . At the few photon level, 
EIT is characterised by the propagation of coupled ex- 
citations of atoms and light called dark-state polaritons 
(DSPs) [2]. These long-lived quasi-particles simultane- 
ously possess the properties of both the photonic and the 
atomic degrees of freedom which can be interchanged in 
a fully coherent and reversible process. As coupled ex- 
citations of light and matter, DSPs acquire many new 
properties from the atomic system. This has been used, 
for instance, to slow down or even stop light [2H6], to im- 
print a magnetic moment onto light fields [7] , or to realize 
giant electro-optical effects in atomic vapours [8]. Qual- 
itatively new many-body effects occur when coupling 
light to strongly-interacting atomic systems. Polariton- 
polariton interactions mediated by the atomic admixture 
can lead to highly nonlinear [9HTT] and nonlocal optical 
effects [12] as well as the emergence of correlations in both 
the atomic and the light fields [I3][l4]- The ability to pro- 
duce and coherently control the propagation of quantum 
fields using interacting dark-state polaritons will open 
up new applications including few photon nonlinear op- 
tics ISIEIHg, non-classical light sources [TQl [T3l HfM9] . 
photonic quantum logic gates [T9H2T] . the generation of 
squeezed states of atomic ensembles [22] and new types 
of strongly-interacting quantum gases [23H25] . 

Here we experimentally study both light and matter 
aspects of strongly interacting DSPs in a dense atomic 
gas involving highly excited Rydberg states. An impor- 
tant feature of this system is that the Rydberg blockade, 
in which a single excitation shifts the Rydberg states of 
nearby atoms out of resonance within a critical distance 
R c [26 , 27 , gives rise to dissipative hard-core interactions 
and ultimately to correlations between DSPs [14]. In con- 
trast to the usual approach, in which EIT is observed op- 




FIG. 1. Experimental setup, a, Level scheme used to 
observe electromagnetically induced transparency in a Ryd- 
berg interacting gas. b, Cross-section of the atomic cloud 
(blue) showing the effect of Rydberg blockade on polariton 
statistics and nonlinear absorption for small and large po- 
lariton densities, c, Schematic of the experimental setup. 
The atomic cloud is illuminated by the probe field yielding 
an absorption image with a local transparency due to the 
coupling beam. Electrodes (segmented rings) and a micro- 
channel-plate detector allow for field ionisation and detection 
of Rydberg atoms. 



tically, we explore a new approach which benefits from a 
high detection sensitivity to Rydberg atoms and directly 
probes the matter-part of the polariton wavefunction. In 
this way we extend upon recent work which has shown 



the development of nonclassical photon correlations in 
the transmitted light field [lOj [13l [19] . Since the map- 
ping between the properties of the light field and the 
atomic system is in general non-trivial [10], it is impor- 
tant to study both aspects to obtain a complete picture 
of EIT in strongly interacting systems. Our experiments 
cover four orders of magnitude in atomic density which 
we use to observe the nonlinear optical response and tem- 
poral dynamics of polar itons, as well as to investigate the 
transition to sub-Poissonian polariton statistics when en- 
tering the strongly interacting regime. Using a state-of- 
the-art theoretical approach for the coupled atom-light 
system we find good agreement with all aspects of the 
experiment except for the level of sub-Poissonian fluctu- 
ations. We observe atomic correlations that are signifi- 
cantly stronger than can be accounted for by assuming 
the propagation of classical light fields alone. This indi- 
cates that the recently observed modification of photon 
statistics of propagating light fields due to Rydberg in- 
teractions [10] can cause a significant back-action on the 
Rydberg atom statistics. 

The basic principle of EIT involves a three-level sys- 
tem (Fig. [T^l) : two long-lived states \g) and \R) are 
coupled via a short-lived state \e) which spontaneously 
decays with rate Y. Strong laser coupling Vt of the 
\e) — >• \R) transition creates an Autler-Townes doublet 
of dressed states, such that the transition amplitudes for 
the \g) — > \e) resonance interfere destructively, giving rise 
to a narrow transparency resonance. This is accompanied 
by the evolution of the atomic system into a dark-state: 
a stable superposition of \g) and \R), which is decoupled 
from the probe field and is thus immune to spontaneous 
decay. In the weak probe limit, a single probe photon 
will be partially converted into a collective excitation of 
the medium to become a dark-state polariton [2], 
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where the collective quantum states involving N atoms 
and the light field are \g^) = \ gi ...g N )\l) and \RW) = 

l/V^VEili \gi.-Ri-.gN)\0). The mixing ratio tan(0) = 
yj pacT I ft determines the properties of the polariton de- 
pending on the atom density p, the resonant scattering 
cross-section a = 3A 2 /27r, and the speed of light c. The 
DSP propagates through the atomic medium with negli- 
gible loss and reduced group velocity v g ~ cj tan 2 (#). 

To understand how interactions between DSPs arise we 
provide the following simple physical picture (Fig. [TJd). 
For weak probe intensities and for low atomic densities 
each incident photon is converted into a DSP which prop- 
agates freely through the EIT medium before being co- 
herently converted back to the original optical mode. If 
the light field is subsequently detected on a CCD cam- 
era this is observed as transparency of the atomic cloud. 
For increasing atomic density the group velocity drops 
(6 — >• 7r/2), effectively compressing the probe field inside 
the atomic cloud and increasing the density of DSPs p DSP 
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FIG. 2. Nonlinear optical response of the Rydberg EIT 
medium. Absorption images of the atomic cloud for three 
expansion times with a probe intensity of I p — 5 /iW/cm 2 
and peak densities corresponding to a, 5 x 10 9 cm -3 , b, 
5 x 10 10 cm" 3 and c, 3 x 10 12 cm" 3 . The EIT region il- 
luminated by the coupling laser is indicated by white ellipses, 
d Measured absorption scaled to the two level response as a 
function of peak atomic density. The dotted horizontal line 
shows the expected low density EIT absorption due to the 
finite laser linewidths. The solid and dashed curves are the 
results of the Monte-Carlo simulations of the coupled atom- 
light system with ^y gR =1.7 MHz (solid) and j g R=2.6 MHz 
(dashed). 



until several begin to occupy the same volume. Since a 
single polariton includes a large Rydberg state compo- 
nent, it interacts strongly with the Rydberg states of 
nearby atoms which become shifted out of resonance, 
breaking the EIT condition for additional excitations and 
leading to strong inelastic scattering of excess photons 
from the cloud [13] . This imposes a critical density of 
DSPs p c « (47ri?;?/3) _1 corresponding to approximately 
one per blockade volume, beyond which the density of 
DSPs saturates (FigjlJ) lower panel). 

Our experimental setup is shown schematically in 
Fig. lb. We use 87 Rb atoms with states \g) = |55 , i/ 2 , F = 
2), fe) = |5P 3/2 ,F = 3) with Y = 6.1 MHz, and the 
Rydberg state \R) = |55Si/2) which has a lifetime of 
« 80 ps. We prepare a cigar shaped cloud of approxi- 
mately 5 x 10 4 atoms at a temperature of 5 pK in an 
optical trap. The atomic cloud is uniformly illuminated 
from the side by a circularly-polarized probe beam with 
intensity I p = 5 pW/cm 2 and the transmitted light is 
imaged onto a CCD camera with an optical resolution 
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of « 12 /im (Rayleigh criterion). The coupling laser is 
counter aligned to the probe beam and passed through 
a diffractive optical element to create an approximately 
uniform intensity profile with Q = 5.1 MHz over an el- 
liptical region of ~ 65 x 130 jam 2 . The frequencies of 
both lasers are tuned to the atomic resonances. The 
measured linewidths (including both decay and dephas- 
ing terms) for the probe transition and the two-photon 
transition are j eg ~ 6.4 MHz and j 9 r ~ 1.7 MHz 
respectively, which gives an EIT coupling parameter 
C = Q 2 1 '{pfeglgRj = 2.4 and an EIT resonance width 
of w = (1 + C)j gR = 5.8 MHz. The anticipated block- 
ade radius is R c — (2Ce/w) 1 ^ 6 « 5 /im with the van-der- 
Waals interaction strength Cq « 50 GHz /im 6 . To control 
the atomic density before pulsing on the probe light we 
switch off the optical trap and vary the time-of-flight be- 
tween 20 /is and 4.5 ms, varying the peak density (along 
the line-of-sight) between p « 10 13 — 10 9 cm -3 . This is 
coupled to a change in the cloud length in the probe di- 
rection between w 2 /im and 100 pm and a change in the 
optical density of the cloud between 3 (resolution limited) 
and 0.1. 

Figure [2] shows absorption images of the atomic cloud 
for three expansion times corresponding to three differ- 
ent peak atomic densities. For long time-of- flights (cor- 
responding to low densities p < 10 10 cm -3 ) the cloud 
has expanded considerably, such that the excitation vol- 
ume is defined primarily by the intensity profile of the 
coupling laser. Here we observe a high degree of trans- 
parency, which can be attributed to the low density of 
DSPs. By analysing the images we measure the absorp- 
tion in both the EIT region and the background region 
corresponding to the two-level optical response (with sus- 
ceptibility X2l)- At low densities we find the optical re- 
sponse is linear corresponding to a scaled optical suscep- 
tibility X0/X2L = 0.3 ±0.1 which is consistent with the 
expected on-resonant optical susceptibility in the non- 
interacting EIT regime X0/X2L — 1/(1 + C) « 0.29. At 
densities above p ~ 1 x 10 10 cm -3 we start to observe a 
decrease in the transparency due to interactions. Above 
p « 4 x 10 11 cm -3 the EIT spot vanishes almost com- 
pletely. 

The full nonlinear optical response as a function of 
peak atomic density is shown in Fig. [2]i. We observe 
a transition from transparent to fully absorbing by in- 
creasing p. This characteristic behaviour resembles that 
obtained by full many-body calculations of the optical 
susceptibility in a Rydberg blockaded gas [28]. To make 
a quantitative comparison between theory and experi- 
ment we use an optimised version [29 of the rate equa- 
tion model originally introduced by Ates et al. [30 . We 
assume random atom positions distributed according to 
the measured Gaussian cloud shape and use Monte-Carlo 
sampling to solve for the steady state populations in 
states \e) and \R) for each atom, taking into account 
the level shifts due to interactions with nearby Rydberg 



atoms. Comparing the simulations to the experiment 
we find that the effect of nonlinear absorption as the 
probe propagates through the cloud must be accounted 
for. Therefore, to solve for the steady state populations 
for any given atom, we first solve for the local probe in- 
tensity which is attenuated by all preceding atoms which 
in turn depends on their steady state populations. The 
coupled local probe intensity and steady state popula- 
tions are then determined recursively. 

We find good agreement between the measured optical 
response and the results of the Monte-Carlo simulations 
using the calculated interaction strength and all other in- 
put parameters determined from independent measure- 
ments (Fig. [2]i-solid curve). The simulations qualita- 
tively reproduce the observed scaling of the nonlinear 
absorption, however at intermediate and high densities 
the scaled absorption is underestimated slightly. The 
discrepancy could be explained by the combination of 
experimental uncertainties in the measured dephasing 
rates, additional density-dependent dephasing of the ^in- 
coherence or by decay to nearby Rydberg states which are 
no longer directly coupled to the light field. In the low 
density limit the simulations with the measured value of 
jgH = 1.7 (solid line) are consistent with the experimen- 
tal data. For comparison we show a fit to the data assum- 
ing a larger overall dephasing of 7^=2.6 MHz (dashed 
curve) . 

In our experiments the DSPs are almost entirely 
matter-like (in eq. cos 2 (<9) w 10" 3 - 10" 7 ). There- 
fore, measuring the Rydberg population serves as a pro- 
jective measurement of the number of polaritons inside 
the cloud. This is possible by field ionizing the Ry- 
dberg states and subsequently detecting the individual 
ions on a microchannel plate (MCP) detector which pro- 
vides high time resolution and near single particle sen- 
sitivity (Fig. [3^l) . We apply an electric field just above 
ionization threshold to field ionize and guide the resul- 
tant ions to maximize the spread in arrival times. In a 
typical experiment we detect « 10 ions. The probability 
for two ions to arrive simultaneously such that they are 
miscounted as a single event is less than 1%. 

Figure ^jp shows the measured time evolution of the 
number of detected ions during a 2.2 ps probe pulse in 
the non-interacting regime. At short times we observe 
a rapid increase of the detected ions reflecting the en- 
try of DSPs into the atomic cloud. The observed rise 
time of ~ 100 ns is compatible with the spectral width 
of the EIT-resonance. After this time the number of de- 
tected ions approaches a steady state value of around 
(Anions) ~ 10, which, accounting for our ion detection ef- 
ficiency (rj « 0.4) corresponds to « 25 polaritons inside 
the cloud. After the pulse the Rydberg state popula- 
tion decreases rapidly, which is expected as the polaritons 
leave the atomic cloud and are coherently mapped back 
onto the light field. By analysing the images we confirm 
that the cloud remains mostly transparent until the elec- 
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FIG. 3. Polariton dynamics observed during a 2 /xs 
probe pulse measured via the Rydberg state popula- 
tion, a, Typical ion signal showing the detection of individual 
ions (marked by triangles). The dashed horizontal line indi- 
cates the threshold voltage of -0.02 V below which we count 
voltage spikes as ion detection events, b, Measured polariton 
evolution with each data point computed from the average of 
20 runs of the experiment. The vertical error bars represent 
the statistical error of the mean. The solid line is a fit to the 
data using the numerical solution to the optical Bloch equa- 
tions, including a coupling to both Rydberg Zeeman sublevels 
rrij — ±1/2 which accounts for the slow evolution and resid- 
ual Rydberg population after 2 /xs. The dashed red line shows 
the population in the resonantly coupled Zeeman state alone. 



trie field is switched on. Over the full pulse duration we 
estimate a total of ~ 1500 polaritons pass through the 
atomic cloud. After the pulse we observe a residual pop- 
ulation of about 30% of the peak number of ions and a 
slower time evolution. The dashed line in Fig. ^jp shows 
the expected evolution found by solving the three-level 
optical Bloch equations. The additional slow rise during 
the pulse and the longer tail observed in the experiment 
can be attributed to laser coupling to additional Zee- 
man sublevels assuming polarization dependent coupling 
strengths (Fig. [3)3- solid line). 

Now to investigate the appearance of correlations in 
the strongly interacting regime, we map out the distri- 
bution of DSPs inside the cloud as a function of atomic 
density. After 2 ps of probe laser illumination we turn 
on the ionization field and count the number of detected 
ions. This is repeated between 150-300 times for each 
time- of- flight corresponding to different peak densities of 
the cloud. From these data we construct histograms of 



the ion number distribution (Fig. |4pi) which is used to 
extract the mean and variance of the number of polari- 
tons. 

In Fig. [4}d we show the measured mean ion number 
(Anions) which is proportional to the number of DSPs in- 
side the cloud. The low-density behaviour is determined 
by the changing number of atoms in the excitation vol- 
ume as a function of time-of- flight (dashed-line). For 
p > 10 11 cm -3 we observe a rapid decrease of (A^ ons ) 
due to probe laser attenuation and polariton-polariton 
interactions. The solid-line is the results of the Monte- 
Carlo simulation using the same parameters as for Fig. [2] 
To obtain the best agreement we adjust only two param- 
eters in the simulation: the detection efficiency rj and 
the effective area of the coupling laser ellipse. Optimal 
agreement is found for an elliptical coupling region with 
minor and major axes of 65 /im and 130 /im respectively, 
and rj = 0.4 which is consistent with our experimental 
geometry. 

Simultaneous to the onset of the polariton blockade 
we also observe a significant narrowing of the statistical 
distribution of the detected ions (Fig. [4^l). We quantify 
this using the Mandel Q parameter 



Q 



(NLA ~ i^i ) s 

(A" ions) 



(2) 



which compares the width of an observed distribution to 
a Poisson distribution with the same mean. For a Pois- 
sonian process Q = 0, whereas a sub-Poissonian process 
yields Q < reflecting the emergence of spatial and tem- 
poral correlations [3TH33] . 

Figure 4c shows the measured Mandel Q parameter 
as a function of the peak density of the atomic cloud. 
In the low-density limit, polariton-polariton interactions 
can be neglected and the measured ion number distri- 
butions are close to Poissonian with Q « which is ex- 
pected for a coherent input light field. Increasing the 
density we observe a clear transition to sub-Poissonian 
statistics with Q reaching a minimum value of —0.32 with 
a statistical uncertainty of 0.04 calculated in the range 
2 x 10 11 cm -3 < p < 2 x 10 12 cm -3 . Including the effect 
of imperfect ion detection indicates that the true distri- 
bution of DSPs is much narrower corresponding to strong 
spatial and temporal correlations (since Qdsp = Q/v)- 

Despite the good agreement of the Monte- Carlo model 
with both the optical response data (Fig. [2]i) and the 
mean number of polaritons (Fig. [4)3), it clearly fails 
to describe the experimentally measured Q parameters 
which are much lower than predicted (Fig. [4^). We at- 
tribute this to the neglected effect of photon correlations 
(the Monte-Carlo algorithm assumes a classical field at 
each step characterised by a Rabi frequency Q p ) which 
can build up as a consequence of dissipative polariton- 
polariton interactions [TP] . 

A complete and rigorous many-body treatment of this 
problem, including the effects of the quantized electro- 
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FIG. 4. Statistical distribution of the detected ions 
during a 2 ps probe pulse showing the transition to 
sub-Poissonian statistics, a, Typical ion number his- 
togram for p = 5 x 10 11 cm -3 with Q — —0.3 compared 
to a Poisson distribution with the same mean (solid line), b, 
Measured mean number of ions as a function of atomic den- 
sity which decreases above p = 10 11 cm -3 due to blockade. 
The solid line is the result of the Monte-Carlo calculations 
with r] — 0.4 and the dashed line shows the Rydberg number 
neglecting interactions and probe absorption, which is propor- 
tional to the total number of atoms in the excitation volume, 
c, In the high density regime the Mandel Q parameter falls to 
sub-Poissonian values. The solid line is the result of Monte- 
Carlo calculations using the same parameters as in Fig.[2]i and 
Fig. [4}x The dashed curve is derived from the hard sphere 
model assuming an initial Poissonian photon number distri- 
bution with pcrit = 1 x 10 12 cm -3 . The statistical errors are 
estimated by dividing the data for each density into 3 sub- 
sets of 50 experimental runs each and then calculating the 
standard error of the mean Q parameter. 



magnetic field is still lacking. However, we can make a 
simple estimate for the possible effects of photon corre- 
lations by assuming a hard-sphere model for polariton 
interactions. We discretize the excitation volume into el- 
ements corresponding to the size of a blockade sphere and 
assume an initial Poisson distribution for the number of 
photons per element with a mean of (N p h) ~ -^p/p c [13] - 
For Nph > 1 we assume that one DSP is created with 
unity probability and the excess photons couple to bright 



states which do not include a Rydberg state compo- 
nent. Therefore the corresponding DSP statistics range 
from Poissonian-character for (N p h) <C 1 (limited by the 
classical fluctuations of the laser field) to strongly sub- 
Poissonian for (N p h) ^> 1 above the critical density at 
which the DSP density saturates. The atomic density at 
which correlations become important is when (N p h) — 1, 
which for our experimental parameters corresponds to 
Pcrit ~ 1 x 10 12 cm -3 . According to this simple model 
we find Q ~ 7yexp(— p/p cr it) — V> which is in closer agree- 
ment with the data, confirming that the development of 
photon number correlations within the cloud plays an 
important role on the DSP statistics (Fig. [ij:). 

For the largest densities we find that the measured 
Q values increase again to super- Poissonian values. We 
have not yet identified the cause of this increase, but we 
note that it co-incides with the point when the cloud 
becomes quasi- ID with respect to the Rydberg excita- 
tion which might modify the effects of interactions on 
the propagating light field. Alternatively it could be due 
to the production of a small number of prompt ions (be- 
low our detection sensitivity) [34] which can dramatically 
affect the polariton statistics. Spurious ions could also be 
present in recent experiments studying the effects of Ry- 
dberg blockade on photon statistics [lOj [THJ [19] but to 
our knowledge their effects have not been investigated so 
far. 

Our experiments provide a direct observation of 
strongly-interacting DSPs realized using EIT in an ul- 
tracold Rydberg gas. Strong polariton-polariton interac- 
tions are observed through the density dependent opti- 
cal response and the polariton counting statistics, which 
become significantly sub-Poissonian at high densities, re- 
flecting the emergence of spatial and temporal correla- 
tions between polaritons. The method presented here 
in principle allows for studying the full polariton count- 
ing statistics, thereby also giving access to higher order 
correlations [31]. Two-photon correlations have recently 
been observed in the second order intensity correlation 
function of light g^ 2 \r) retrieved from small Rydberg en- 
sembles p~8j [19] , and via light propagation through dense 
atomic clouds [10 . We point out however that in general 
there is no direct mapping between the correlations on 
the light field and the atomic correlations. For instance, 
in a dense but optically thin medium, excess photons may 
not be efficiently scattered out of the forward direction 
thereby washing out the contrast, despite the pres- 
ence of strong atomic correlations. Alternatively, out- 
side of the blockaded regime, spin- wave dephasing due to 
atomic interactions can alter the statistics of the emitted 
light p~8j [35j [36], even though the atomic correlations 
could remain nearly classical. A complete description 
of EIT in strongly-interacting media remains a theoreti- 
cally challenging problem [10] , therefore it is necessary to 
study both aspects and to provide more stringent tests 
to theoretical models. Ultimately this would open an av- 
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enue to on-demand sources of single photons and other 
nonclassical states of light. 

Our experiments can be interpreted as Rydberg dress- 
ing of photons, in which photon-photon interactions are 
mediated by the Rydberg state admixture. By detuning 
from the intermediate state, the dissipative character of 
the polariton-polariton interactions could be suppressed 
while maintaining the huge nonlinearities which would 
give rise to elastic scattering [14]. As the nonlinearity 
can be controlled via the choice of the Rydberg state, 
repulsive, attractive and long-range interactions become 
possible. This could pave the way to create and study 
new types of quantum gases such as dipolar Bose- 
Einstein condensates of dark-state polaritons [241 125], 



METHODS 

Coupled atom-light simulations. To model the ex- 
periments we use an optimised rate equation model [29] . 
We simulate the full number of atoms within the 
coupling volume (N « 5000 — 25000) assuming ran- 
dom positions distributed according to the measured 
Gaussian cloud shape. Monte-Carlo sampling is used 
to solve for the steady state populations for each atom, 
taking into account the level shifts due to interactions 
with nearby Rydberg atoms. The probability for an 
atom to change its state is determined by a stochastic 
process. In each step an atom is picked randomly and 
the steady state populations are calculated, taking into 
account level shifts due to the interaction with nearby 
Rydberg atoms. To include the effects of probe laser 
attenuation we first solve for the local Rabi frequency 
Qp^ which is attenuated by all other atoms j < i in a 
cylindrical volume with cross-section A. The degree of 
attenuation is given by flp +1 ^ = Qp^ exp (— x^) where 
is proportional to the |e) state population for atom 
j. We recursively determine the local Qp^ starting with 
j = 1 which sees the total incident probe Rabi frequency 
ttp 1 ^. We then solve for the steady state populations 
iteratively until j = i. The outcome of this procedure 
does not depend on the exact choice of A (typically 
3 — 10 fim 2 ) as long as A is small compared to the size 
of the atomic cloud. 

Counting polaritons. The number of polaritons inside 
the cloud is measured by field ionizing the Rydberg states 
and guiding the ions to an MCP detector. The field is 
switched above the field ionization threshold within 300 
ns and the resultant ion trace is recorded on a fast os- 
cilloscope. Each ion creates a voltage spike with a width 
of ~ 1.5 ns with arrival time distributed over ~ 600 ns. 
A/i ons is found by fitting each trace with the MCP transfer 
function using sparse spikes deconvolution and counting 



extrema for pulses with amplitudes less than the thresh- 
old voltage of -0.02 V. From the pulse-height-distribution 
we infer a counting efficiency of 0.9, which combined with 
the MCP detection efficiency, angular responsivity and 
guiding efficiency gives an estimated overall efficiency 
which is in good agreement with the best fit to the data 
of r] ~ 0.4. For the low numbers of ions produced in our 
experiment detector saturation can be safely neglected. 
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